Self-consistent description of simplified pseudospin-electron model 
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A method of the self-consistent calculation of the thermodynamical and correlation functions is presented. This 
approach is based on the GRPA (generalized random phase approximation) scheme with the inclusion of the mean 
field corrections. The numerical research shows that interaction between the electron and pseudospin subsystems 
leads to the possibility of the dipole (pseudospin) and charge-density instabilities as well as phase separation into 
the uniform and/or chess-board phases. 



The aim of this article is to calculate the corre- 
lation functions, pseudospin and particle number 
mean values as well as the grand canonical poten- 
tial for the simplified (Z7 = 0) pseudospin-electron 
model (PEM) within the one self-consistent 
approach. 

PEM includes the local interaction of the con- 
ducting electrons with some two level subsys- 
tem described by pseudospins S z = ± 1/2 (e.g an- 
harmonic vibrations of the apex oxygen ions in 
YBaCuO-type crystals) 0. On the other hand, 
PEM can be transform into the binary alloy type 
model as well as into the Falikov-Kimball model 
(FKM) and they can be also studied analytically 
within the framework of the developed scheme. 

The model Hamiltonian is the following: 

H=Y,H i0 +Y J t ij cj lT c ja , (1) 
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where an interaction with pseudospins (g-term) 
placed in some longitudinal field h (chemical po- 
tential for ions in FKM) are included in the single- 
site part; fi is the chemical potential. The sec- 
ond term in the Hamiltonian describes an electron 
hopping from site to site. 

The calculations are performed in the strong 
coupling case (g 3> t) using single-site states 
as the basic one. A formalism of the electron 
annihilation (creation) operators = a a P^, 
a-ia = CiaPf (P[° — \ i <Sf ) acting at a site with 
the certain pseudospin orientation is introduced. 

Expansion of the calculated quantities in terms 



of the electron transfer leads to an infinite series of 
terms containing the averages of the T-products 
of the die hia operators. The evaluation of such 
averages is made using the corresponding Wick's 
theorem. The results are expressed in terms of 
the products of the nonperturbed Green's func- 
tions and averages of the projection operators 
P± which are calculated by means of the semi- 
invariant expansion 0j. 

The calculation of the correlation functions is 
performed within a self-consistent scheme in the 
framework of the generalized random phase ap- 
proximation (U with the inclusion of the mean 
field type contributions coming from the effective 
pseudospin interaction via conducting electrons, 
i.e. all zero-order correlators (second order semi- 
invariants) as well as pseudospin mean value are 
calculated by the mean-field Hamiltonian: 
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Diagram equation on pseudospin correlator 
(S x S z ) g is following: 




This equation differs from the one for the fsing 
model in MFA by the replacement of the exchange 



interaction by the electron loop (which describes 
an interaction between pseudospins via conduct- 
ing electrons). In the analytical form its solution 
is equal 
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Equation for pseudospin mean value in the uni- 
form case ((Sf } = (S z )) is following: 

1 f l+e _/3ei 
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We also calculate particle number mean value and 
grand canonical potential |3j. All quantities can 
be derived from the grand canonical potential by 
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which show the thermodynamical consistence of 
the proposed approximations. 

The analysis of the pseudospin correlator tem- 
perature behaviour shows that high temperature 
phase become unstable with respect to fluctua- 
tions with wave- vector q — (tt, 7r) (for some model 
parameters values) that indicates the possibility 
of the phase transition into a modulated (chess- 
board) phase. On the other hand, the phase 
transition between the uniform phases with dif- 
ferent pseudospin mean values (bistability) is pos- 
sible g. 

From the comparison of the grand canonical 
potential f2 values for uniform and chess-board 
phases (^=const regime), the (/i, h) phase dia- 
gram is obtained (Fig. la). 

Chess-board phase exists as intermediate one 
between the uniform phases with different (S z ) 
and (n) values. The transition between different 
uniform phases (bistability) is of the first order 
(Fig. la, dashed line), while the transition be- 
tween the uniform and modulated ones is of the 
first (dotted line) or second (solid line) order. 

On the other hand, the minimum of the free 
energy F=fl+fj,N is the equilibrium condition in 
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Figure 1. Phase diagrams: (a) // — h, (b) n — h. 
I - uniform phase, II - chess-board phase, PS - 
phase separation area. (T/g — 0.005, t/g = 0.2) 

the n=const regime. In this regime the first order 
phase transition with a jump of the pseudospin 
mean value accompanied by the change of elec- 
tron concentration transforms into a phase sepa- 
ration into the regions with different phases (the 
uniform and chess-board ones) and with differ- 
ent electron concentrations and pseudospin mean 
values (Fig. lb). 

On the basis of the presented scheme the ther- 
modynamics of phase transitions, possibility of 
the phase separations as well as appearance of 
the chess-board phase have been investigated. 
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